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Quasistatic  and  Dynamic  Growth  of  Microscale  Spherical  Voids 


T,  C,  Tszeng 

Berkeley  Materials  Research 
Berkeley,  CA  94704,  USA 


Abstract 

This  study  examines  the  quasistatic  and  dynamic  growth  of  microscale  spherical  voids.  A 
general  dynamic  model  of  void  growth  is  first  developed  by  considering  conservation  of  local 
energy,  including  external  work,  surface  energy,  kinetic  energy,  clastic  strain  energy,  and  plastic 
dissipation  associated  with  growing  void  driven  by  hydrostatic  tensile  stress.  It  properly  accounts 
for  material  compressibility  and  limited  size  of  elastic  deformation  in  dynamic  growth  of  spherical 
voids.  A  closed  form  of  pressure-void  size  relation  is  obtained  for  quasistatic  growth  of  spherical 
voids  embedded  in  non- strain  hardening  materials.  Critical  condition  is  identified  for  unbounded 
growth,  which  depends  strongly  on  initial  void  size.  For  suberitical  loading,  a  void  can  only  grow  to 
a  very  limited  size.  An  interesting  crossover  phenomenon  is  identified  pertaining  to  the  influences 
of  yield  stress  on  void  growth  rate. 

L  Introduction 

The  phenomena  of  nudeation,  growth  and  coalescence  of  voids  arc  the  fundamental 
mechanism  in  ductile  fracture,  Void  nudeation  can  occur  through  decohesion  between  a  ductile 
matrix  and  second-phase  particles  (Needleman,  1987),  or  fracture  of  brittle  particles  in  a  ductile 
matrix.  There  arc  indications  that  lattice  defects  in  the  length  scale  of  nanometer  can  also  be  the 
embryo  of  void  nudeation  (Lubarda  et  al.,  2004).  Once  voids  have  been  nucleated,  they  grow 
through  localized  plastic  deformation.  There  are  two  broad  subjects  pertaining  to  void  growth: 
threshold  condition  and  dynamics  of  void  growth.  On  the  first  subject,  many  arc  also  focused  on 
bifurcation  and  instability,  Needleman  (1972)  established  a  model  to  describe  the  growth  of 
spherical  void  in  elastic-plastic  medium.  Chung  et  al.  (1987)  investigated  cavitation  bifurcation 
under  a  remote  hydrostatic  tensile  loading,  and  concluded  that  cavitation  pressure  was  of  the  order 
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of  Young's  modulus.  Huang,  Hutchinson,  and  Needleman  (1991)  addressed  quasistatic  cavitation 
instabilities  in  elastic-plastic  materials  subjected  to  niultiaxial  axi symmetric  stressing  and  found 
that  the  condition  for  unstable  void  growth  depends  on  the  attainment  of  a  critical  value  of  the  mean 
stress.  Tvergaard  and  Hutchinson  (1993)  extended  this  analysis  to  a  power  hardening  elastic-plastic 
material  and  also  concluded  that  the  occurrence  of  cavitation  instability  depends  primarily  on  the 
level  of  the  mean  stress.  Hou  ct  a 1.  (1999)  considered  effects  of  surface  energy  on  quasistatic 
growth  of  spherical  and  cylindrical  voids,  Tvergaard  and  Hutchinson  (2002)  found  that  initial  void 
shape  had  little  effect  on  the  critical  stress.  Thus  it  appears  that  both  deviations  from  the  hydrostatic 
state  and  deviations  from  the  spherical  initial  void  can  be  neglected  in  the  study  of  void  growth  (Wu 
et  a  l,,  2003), 

Hunter  and  Crozicr  (1968)  and  Crozicr  and  Hunter  (1 970)  addressed  the  dynamics  of  void 
growth  by  considering  compressibility  in  elastic-plastic  material,,  and  obtained  similarity  solution 
for  steady-state  growth.  Most  published  studies  on  the  dynamics  of  void  expansion  take  advantage 
of  material  incompressibility,  and  primarily  on  spherical  voids.  Rice  and  Tracey  (1969)  employed 
Rayleigh-Ritz  approximation  to  a  variational  principle  for  the  flow  field  in  an  elastically  rigid  and 
incompressible  plastic  material  containing  a  void,  subjected  to  a  remotely  uniform  stress  and  strain 
rate  field.  For  spherical  voids,  Carroll  and  Holt  (1972)  addressed  the  dynamic  collapse  of  a  void  in 
an  ideally  plastic  material  by  using  a  hollow  sphere  model.  Johnson  (1981 )  extended  this  approach 
to  rate-dependent  materials  under  hydrostatic  tensile  loading.  The  dynamic  void  collapse  and  void 
growth  in  single  crystals  under  uniform  farfield  stresses  were  studied  analytically  by  Memat-Nasser 
and  Hori  (1 987),  Ortiz  and  Molinari  (1992)  considered  the  dynamic  growth  of  a  single  void  in  a 
power-law  hardening  material  from  an  energetic  viewpoint.  Cortes  (1992a)  extended  the  spherical 
shell  approach  to  examine  void  growth  under  combined  hydrostatic  and  dcviatoric  stresses  in  rigid 
plastic  materials;  Cortes  (1992b)  addressed  strain  hardening,  rate  sensitivity  and  thermal  effects  of 
the  Johnson-Cook  type,  but  elasticity  was  neglected  because  of  the  assumption  that  plasticity  would 
quickly  dominate  the  volume.  Wang  (1994)  further  extended  this  approach  to  model  void  growth  in 
rate-dependent  materials.  Wu,  Ramesh  and  Wright  (2003)  included  the  influence  of  thermal 
diffusion  and  thermal  softening. 


There  are  also  studies  based  on  micromechanics  of  dislocations  pertaining  to  stability  and 
dynamics  of  void  growth  (e.g.,  Stevens  el  al.,  1972;  Kameda,  1989;  Lubarda  et  at.,  2004;  Ohashi, 
2005;  Ahn  et  al.,  2006,  2007).  These  studies  are  mostly  focused  on  high  rate,  high  magnitude 
loading  {e.g.,  shock)  in  which  deformation  mechanisms  associated  with  void  formation  maybe 
different  from  regular  loading  at  low  rate.  Based  upon  dislocation  micromcchanics,  the  present 
author  has  studied  the  threshold  condition  in  general  macroscopic  stress  state  (Tszeng,  2007a),  and 
dynamic  model  of  transient  growth  (Tszeng  2007b)  associated  with  cylindrical  voids.  Significant 
progress  has  also  been  achieved  recently  on  void  nucleation,  growth  and  coalescence  by  the 
technique  of  atomistic  and  Molecular  Dynamics  simulations  (e.g.,  Rudd  and  Belak,  2002;  Scppala 
et  al,  2005). 

The  subject  of  dynamic  expansion  of  cylindrical  void  has  been  mostly  avoided  in  the  past. 
By  assuming  incompressibility,  it  can  be  shown  that  total  kinetic  energy  for  an  expanding  spherical 
void  is  K  =  2xpR\>2  (Ortiz  and  Molinarf  1992;  Wu  et  al,  2002),  where  R  is  void  radius  and  v  = 
dR/dt  is  the  growth  rale.  If  the  same  procedure  is  applied  to  cylindrical  void,  the  kinetic  energy 

R 

would  be  K  =  xpR2v2  In  —  „  where  is  the  outer  radius  of  deformation  zone  (Tszeng,  2007c), 

R 

Incompressibility  requires  R.,  ->  and  therefore  kinetic  energy  K  is  infinite.  In  another  way,  stress 
equilibrium  in  radial  direction  of  a  cylindrical  void  can  be  written  as  (Crozier  and  Hunter,  1970; 
Tszeng,  2007c,  d); 

d2R  .  dR 

dr  +  dj  *  (!) 

dr  r  r 

where  et,  is  radial  stress,  R0  initial  void  radius,  G  shear  modulus,  and  p  mass  density. 

Logarithmic  singularity  results  from  integrating  the  terms  on  the  right-hand  side  of  Hq.  {1 )  to 
infinity.  Its  implication  is  a  zero  growth  rate  for  cylindrical  void.  Obviously,  infinite  deformation 
zone  is  the  underlining  issue  that  prevents  the  nonsingular  solution.  Therefore,  Crozier  and  Hunter 
(1970)  and  Warren  (1999)  declared  that  no  physically  realistic  solution  is  possible  for  steady- 
dynamic  expansion  of  long  cylindrical  cavities  (plane-strain  assumption)  in  an  incompressible 
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elastoplastic  medium,  due  to  logarithmic  stress  divergence  in  the  remote  elastic  field.  Warren 
(1999)  examined  the  effects  of  strain  rate  on  the  model  of  Crozier  and  Hunter  (1970)  for  cylindrical 
void  in  steady  expansion.  Just  recently,  Masri  and  Durban  (2006)  challenged  the  declaration,  of 
Crozier  and  Hunter  (1970),  and  obtained  self-similar  solutions  to  the  model  of  steady-state  growth 
using  hypoelasticity  h  theory.  A  disposable  constant  was  used  to  remove  stress  singularity. 
Unfortunately,  as  pointed  out.  in  Tszcng  (2007c),  there  is  a  severe  internal  inconsistency  that  defeats 
the  validity  of  that  model.  That  is,  stress  singularity  does  exist,  and  the  solution  is  approximate,  at 
the  best.  In  the  case  of  spherical  void,  stress  equilibrium  in  radial  direction  is  (Wu  ct  ah  2003) 


da 


dr 


2  a  RR  +  2RR2  ^ 

—  =  P( - i - 2 

r  r 


R*R2 

5  J 


(2) 


Integrating  the  terms  on  the  right-hand  side  of  Eq.  (2)  to  infinity  docs  not  lead  to  singularity.  The 
nonsingular  solution  is,  however,  an  approximation  valid  for  incompressible  materials.  It  is  believed 
that  considerations  should  be  given  to  limited  deformation  zone  to  circumvent  the  same  issue  that 
prevents  proper  description  of  the  dynamics  of  expanding  cylindrical  voids. 

This  paper  examines  the  quasistatic  and  dynamic  growth  of  spherical  voids  in  sub-micron 
length  scales.  A  general  dynamic  model  of  void  growth  is  first  developed  by  considering 
conservation  of  local  energy,  including  external  work,  surface  energy,  kinetic  energy,  elastic  strain 
energy,  and  plastic  dissipation  associated  with  growing  void  driven  by  hydrostatic  tensile  stress. 
Special  attention  is  placed  on  the  development  of  a  proper  expression  of  kinetic  energy  in  the 
material  surrounding  an  expanding  void.  In  Section  2,  the  governing  equation  of  void  growth  is 
established.  Section  3  examines  quasistatic  growth  and  threshold  conditions.  Lastly,  dynamic 
growth  is  discussed  in  Section  4,  with  special  emphasis  on  steady-state  growth, 

2.  Dynamic  Model  of  Void  Growth 

A  common  difficulty  to  most  dynamic  models  of  void  growth  in  the  past  stems  from  the  fact 
tli at  deformation  extends  to  infinity  for  incompressible  material.  Recognizing  this  fact.  Hunter  and 
Crozier  (1968)  considered  compressibility  in  material  surrounding  an  expanding  void.  They 
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obtained  similarity  solutions  to  the  problem  of  steady  growth  of  spherical  voids.  Obviously,  the 
solutions  are  limited  in  applicability. 

The  present  approach  considers  a  system  that  contains  a  single  spherical  void  embedded  in 
an  infinite  elastic-plastic  matrix.  The  space  is  divided  into  several  zones,  including  void  itself 
incompressible  plastic  zone5  compressible  elastic  zone,  and  undisturbed  material  (Fig.  1).  In  the 
following,  equations  are  formulated  for  elastic  and  plastic  zones  in  the  framework  of  continuum. 
Solutions  are  sought  for  plastic  and  clastic  zones,  respectively,  with  proper  continuity  imposed  at 
the  interface.  In  particular,  the  radial  stress  on  plastic-elastic  interface  is  determined  by  the  solution 
of  elastic  zone.  As  will  be  further  elaborated  in  Section  3,  solution  of  radial  stress  on  plastic-elastic 
interface  as  obtained  by  Hunter  and  Crozicr  (1968)  does  not  reconcile  with  the  solutions  based  upon 
near-incompressible  materials  or  low-rate  growth.  On  the  other  hand,  an  approximate  solution  is 
obtained,  which  is  applicable  for  general  dynamic,  transient  growth. 

2,1  Solution  for  the  plastic  zone 

A  void  has  radius  R  at  current  time  its  surface  moves  outwards  at  radial  velocity  v  (Fig. 

I).  Material  gains  kinetic  energy  associated  with  moving  void  surface.  There  is  energy  dissipation 
due  to  plastic  deformation  and  surface  energy  expenditure  for  expanding  void  surface.  The  law  of 
energy  conservation  states  that. 


W  =  t  +  K  +  ct> , 


(3) 


where  W  is  input  power  done  by  externally  applied  stress/s  train,  f  surface  energy,  K  kinetic 
energy,  and  <h  plastic  dissipation.  Elastic  strain  energy  is  neglected  in  plastic  zone  (will  be 
considered  in  the  outer  elastic  zone;  see  Sections  2.2.  and  2.3).  Eq.  (3)  says  that  surface  energy  and 
kinetic  energy  increase  at  the  expense  of  externally  supplied  energy,  with  additional  energy 
dissipating  in  plastic  deformation. 

The  plastically  deforming  matrix  has  a  radius  D  >  R.  The  incremental  displacement  at.  The 
outer  boundary  of  matrix  material  at  radius  D  in  a  time  interval  dr  is  du,  Assume  the  radial  stress 
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appears  at  r  —  D  is  <rfJ,  the  work  done  in  plastic  zone  is  dW  =  4?iD2crDdu  ,  The  explicit  expression 
of  &D  is  determined  in  Sections  2,2  and  2,3,  By  considering  incompressibility  of  material  between 
R  and  Ds  radial  velocity  at  r  e  (R,  D)  is  (Carroll  and  Holt,  1972;  Ortiz  and  Molinari,  1992) 


V  =(Rirfv, 


(4) 


where  v  =  dRi  dt  is  void  growth  rate,  Radial  velocity  at  the  outer  surface  of  matrix  is  vD  - 


Therefore,  rate  of  energy  input  to  the  plastic  zone  is 


W  =  4xR2<tdv. 


(5.) 


Change  of  surface  energy  is  written  as  dT  =  ydS  where  y  is  surface  energy  density.  Increase  of 
void  surface  dS  =  4 tz(R  +  dR'y  —4 7tR~  =  SxRdR.  Hence, 

T  =  RiryRv  .  (6) 


Material  in  the  vicinity  of  void  surface  gains  velocity  which  is  in  the  form  of  Eq.  (4).  Kinetic  energy 
corresponding  to  an  incremental  mass  at  radius  r  is  dK  =  —  v~dm  =  2h pt2v2.dr  ,  where  p  is  mass 


density.  Total  kinetic  energy  in  plastic  zone  is 


4  -  v" 

K  =  {-xpR’)(— )a, 
3  2 


(7) 


where  a  —  3(1  -  3  //3),  and  ft  =  Di  R .  In  contrast,  the  counterpart  in  a  cylindrical  void  is 

y 

-  V 

K  =  (yr/j/C)(— )a  where  a  =  2  In/?  (Tszeng,  2007c).  In  both  cases,  the  expression  preceding  a 
represents  the  mass  of  void.  Again,  the  geometrical  factor  a  defines  the  influence  of  size  of  plastic 
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zone  on  the  dynamics  of  void  growth.  The  specific  expression  of  D  will  be  discussed  later  in  this 
section.  Change  in  kinetic  energy  between  t-  and  l,  +dt,  is 


dK  =  —  7rpa{R  +  dRy  (v  +  dv'y  —  —  irpaR'v2 . 
3  3 


After  removing  higher  order  terms,  one  obtains  the  expression  of  changing  rate  of  kinetic  energy: 


K  =  2  xpccR 2  v  ■'  +  —  jipaR 3  v — 
3  dt 


(8) 


The  last  term  in  Eq,  (3)  is  plastic  dissipation.  O.  Specific  plastic  dissipation  rate,  W  p  (dissipation 
rate  per  unit  volume),  is 


W  P  =  <76 


(9) 


where  cr  and  &■  arc  the  effective  stress  and  effective  strain  rate,  respectively,  Tn  plastic  zone,  elastic 
strain  is  neglected;  and  total  strain  is  equal  to  plastic  strain;  i.e.,  e  =  sp .  For  simple  strain  hardening 
materials. 


(7  =  <7  yf(s  ") 


(10) 


In  the  present  case  of  spherical  deformation  (Hou  et  ak,  1999), 


£ 


-sr ,  and  sr  =  -2s e  -  \  tn(l  -  — — ; 

3  r' 


R: 


(1 1) 


and 


t'=-i:=2^=2{*Yv- 
r  r  r 


(12) 
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By  Eqs,  (9),  (1 0),  and  (12),  specific  plastic  dissipation  rate  per  unit  volume  is 


^=2<rr/(^)AJ-. 

r  r 


(13) 


Total  plastic  dissipation  rate  associated  with  an  expanding  spherical  void  can  now  be  written  as 


T>  =  f  &Yf{s  )(— ) 2  —  dV  =  Vrv  f  f  (ep)— 

J  r  y  J?  1  p 


(14) 


Hence,  incremental  plastic  dissipation  is  found  to  be 


<&  =  4xR^&¥v, 


(15) 


where  equivalent  flow  stress*  dy ,  is  defined  by 


(iQ 


Eq,  (16)  can  be  further  converted  in  a  form  found  in  literature.  By  using  Eq,  (11)  and  some 

dr  1  ds 


algebraic  operations*  it  can  be  shown  that  —  =  — 

r  2  x  3  _ 


,  Therefore, 


-M 


■i  f{s;:')de 
exp(|f)-l 


exp(-£)-l 


,  where  e.  and  e,.  are  the  effective  strains  at  r  =  R  and  D*  respectively. 


By  Eq,  (11)  again.  £,  =  -2  ln(32rj  /  J?)  =  21n(f? I Ro)  is  the  strain  on  void  surface,  and 
2  R*  —  i?5 

$2  =gy  =—  ln(l  — ^  0 )  c  orresponds  to  the  strain  on  the  boundary  of  plastic  zone.  In  this 

expression,  sy  is  the  strain  at  yielding  (  sy  =  aY  /  E  ).  After  swapping  the  integration  limits*  the 
expression  of  equivalent  flow  stress  in  Eq.  (16)  becomes, 
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f 


|C  Lnl  .'r W .  >  f(  £)  ds 


(17) 


Eqs.  (5),  (6),  (8),  and  (15)  give  the  formula  of  each  term  involving  local  energy 
conservation.  Placing  them  in  Eq,  (3),  one  obtains  the  ordinary  differential  equation: 


(18) 


This  is  the  equation  governing  dynamic  void  growth  of  a  single  spherical  void  embedded  in  infinite 
matrix.  The  radial  stress  on  yield  surface,  ,  is  yet  to  be  identified  in  Sections  2.2  and  2.3. 

In  Eq,  (18),  the  factor  a  is  not  determined  either.  The  boundary  of  plastic  zone  is  defined 
by  the  location  at  which  e  =  ei .  According  to  Eq.  (11),  the  boundary  of  plastic  zone  is  located  at 

radius  D: 


ft  =  ~  =  [1  -  exp(“£j,)l  m t  for  R  »  R0 . 
R  2 


(19) 


3  3 

By  the  approximation  1  -  exp( — eY)  ,  one  has  the  expression  (Huang  el  al.,  1991): 


(20a) 


for  R  »  jV 


As  a  side  note,  Eq.  (20a)  is  actually  valid  for  incompressible  materials;  the  general  expression  ought 
to  be  (Hill,  1989). 


p=[3(i-v)6Yru\ 


(20b) 
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This  ratio  is  presumably  a  material  constant.  Hence. 


«=  3[l-<3(I-v)er)1/3]  (21) 

The  equivalent  flow  stress  &y  of  Eq.  (17)  is  usually  integrated  numerically  providing  the 
function/is  specified.  In  case  a  time’ stepping  method  is  used  for  calculating  the  progress  of  void 
size,  <7y  needs  to  be  updated  accordingly  for  each  lime  step.  As  will  be  shown  later,  radial  stress, 
,  is  also  dependent  upon  current  solution,  and  needs  to  be  updated. 

2,2  Solution  for  the  elastic  zone  -  Steady  state 


The  objective  of  this  section  is  to  determine  the  radial  stress*  crDf  on  yield  surface  at  r  -  D 
in  steady  state  growth.  The  solution  was  obtained  previously  by  Hunter  and  Crozier  (1968).  We 
furnish  the  following  formulation  with  additional  purpose  of  laying  the  ground  for  developing  a 
competing  model  to  be  used  for  general,  transient  void  growth  (Section  2.3).  The  solution  is  first 
sought  for  the  problem  corresponding  to  a  boundary  problem  in  which  constant  velocity  is  specified 
on  the  inner  boundary  of  elastic  zone.  Superposition  principle  is  then  employed  to  determine  the 
radial  stress,  <rD .  Identical  results  arc  obtained  in  an  alternative  approach  by  solving  the  complete 
equilibrium  equation  (Tszeng*  2007 d). 

The  outer  boundary  of  elastic  zone  is  determined  by  the  travel  distance  of  dilatational  elastic 
wave  whose  speed  is 


Ce  ~  ( 


X  +  2G}\ 


(22) 


where  X  is  Lame  elastic  constant,  G  shear  modulus,  and  p  mass  density,  ft  is  assumed  that  void 
grows  from  a  negligibly  small  size.  At  time  f*  the  location  of  outer  boundary  is  given  by  C  =  cj ; 

elastic  zone  is  defined  between  D  and  C  (Fig.  1 ).  Similarity  solution  of  particle  displacement*  w,  in 
the  compressible  elastic  zone  is  obtained  by  Hunter  and  Crazier  (1968)  for  steady-state  growth: 
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(23) 


,  ,2  ,  1  _■), 

u  =  At(-rj-l +-rj  -), 


where  the  dimensionless  similarity  variable  rj  =  —  f  and  A  is  a  coefficient  independent  of  rj. 

CJ 

Particle  velocity  at  a  location  r  e  (D,  C)  in  elastic  zone  is  vp  =  duidt .  By  using  Eq.  (23),  it  can  be 
shown  that  velocity  in  elastic  zone  is: 


vr  =  ^(-l  +  r;-2)  =  A 


„2 


(24) 


Unlike  Hunter  and  Crozier  (1968),  the  present  study  assumes  an  incompressible  plastic  zone, 

Rz 

Particle  velocity  at  the  interface  of  elastic  and  plastic  zone  (r  =  D)  is  =  — -  v ■  The  coefficient  A 

can  then  be  determined  by  imposing  velocity  continuity  at  elastic-plastic  interface;  the  result  is 


A  = 


i  -Fn\ 


(25) 


where  tjR  —  - —  .Therefore,  velocity  field  in  elastic  zone  is 
cJ 


\-rj2  R: 

t2 


On  yield  surface,  von  Mises  yield  criterion  is 


~  r 


By  Hookers  law. 


(26) 


(27) 
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(28) 


E  fv  v  u 

= - (— <t_  +  )  = - tr+— 

"  1-v  E  '  1-v  D 


Replace  Eq.  (28)  in  Eq.  (27),  the  radial  stress  is  oq  written  as 


E  u  l-v 

o'  = - & 

'  1-2 vD  \-2v 


(29) 


By  Eq.  (23)  in  conjunction  with  Eq.  (25),  it  can  be  shown  that 


u 

~D 


1  + v 


^VD 


_2 _ 

fin  r 


(30) 


By  combining  Eqs,  (29)  and  (30),  the  radial  stress  induced  by  a  specified  velocity  on  the  inner 
boundary  of  elastic  zone  is 


=-f^y+ 


1  +  v  P'-jjj 
l-2v  l  +  fii]R 


(31) 


It  can  be  shown  that  radial  stress  of  Eq.  (31)  is  identical  to  that  of  Hunter  and  Crazier  (1968)  in 
steady-state  growth.  For  a  material  subjected  to  macroscopic  hydrostatic  stress  P,  tine  solution  can 
be  obtained  by  simple  superposition;  i.c.,  &D  =P-&r .  Therefore,  the  radial  stress  on  yield  surface, 


fT 


D 


P.l„.  2fffl 

3  '  3  1  —  2v  1  +  flijR 


(32) 


The  expression  of  Eq.  (32)  can  also  be  derived  by  solving  the  stress  equilibrium  in  radial  direction, 
as  shown  m  a  separate  paper  (Tszeng,  2007d), 

By  using  Eqs.  (20b)  and  (22),  Eq.  (32)  is  further  reduced  to 
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(33) 


£>  2  ^2 
etd  =  P  —  ar - <xrv 


3  1  2 


where 


o)  = 


/?(!  +  ffy)  9  I  -  v 


(34) 


The  radial  stress  on  yield  surface,  a D ,  of  Eq,  (33)  is  to  be  used  in  the  governing  equation  of  Eq, 

(18)  for  steady-state  growth.  The  solution  represented  by  Eq.  (33)  is  exact,  subjected  to  the 
assumptions  herein  (Hunter  and  Crazier,  1968).  However,  as  will  be  further  elaborated  in  Section  3, 
it  does  not  reconcile  with  the  solutions  based  upon  incompressible  materials. 

2.3  Solution  for  the  elastic  zone  -  Transient 

The  solution  presented  in  Section  2,2  is  valid  for  steady-state  growth.  A  competing  but  less 
exact  solution  is  now  developed  for  general  dynamic,  transient  growth.  In  the  present  formulation, 
plastic  zone  adjacent  to  elastic  zone  is  assumed  incompressible.  Once  the  macroscopic  stress 
exceeds  a  threshold  condition,  plastic  zone  is  established  immediately.  Compressibility  in  elastic- 
zone  is  considered  in  addressing  the  kinetic  energy  only;  otherwise  incompressibility  is  still 
assumed  in  calculating  elastic  strain  energy.  The  rationale  is  that,  dynamics  of  an  expanding  void 
have  the  greatest  influence  on  kinetic  energy,  and  less  on  elastic  strain  energy. 

The  formulation  is  in  analogy  to  that,  of  Section  2.1,  Local  energy  conservation  is  in  a  form 
similar  to  that  of  Eq,  (2) 


W  =  K  +  '¥ 


(35) 


where  W  is  input  power  done  by  externally  applied  stress/ strain,  K  kinetic  energy,  and  'E  strain 
energy.  Each  term  is  formulated  in  the  following,  Total  energy  input  to  clastic  zone  is 
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W  =  4xR2(P-ctd)v 


(36) 


Total  kinetic  energy  in  clastic  zone  is  K  =  J2jt -pr'v'dr ,  Velocity  field  in  elastic  zone  is  given  by 

Hq.  (26)  for  steady-state  growth.  Nonetheless,  this  equation  describes  a  permissible  velocity  that 
satisfies  boundary  conditions.  As  an  approximation,  we  use  the  same  expression  to  describe  the 
snap  shot  of  velocity  field  in  elastic  zone  at  any  tune  f;  that  is, 


v.  = 


{cjf-r1  R 
{cjf-D2  r2 


■v , 


(37) 


It  is  remarked  that  the  velocity  v  needs  not  to  be  a  constant  anymore.  It  can  be  proven  that  the 
condition  of  velocity  continuity  at  elastic-plastic  boundary  also  satisfies  the  yielding  condition. 
Total  kinetic  energy  is  obtained  for  clastic  zone  between  D  and  cj  .  After  some  algebra  operations, 
the  res li It  is 


where 


(38) 


(3  +  — ffi) 

m+m2 


(39) 


and  normalized  average  growth  rate  v  =  —  \vdi= — .  Actually,  v  =  t]R ,  but  it  is  more  convenient 

cj  J  cj 

to  use  v  for  the  cases  of  transient  growth  being  dealt  with  in  this  section.  Since  Eq.  (38)  is  in  the 
same  form  of  Eq.  (7),  it  is  easy  to  show  that  changing  rate  in  kinetic  energy  is 
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(40) 


K  =  2xpxR V 


4 

+  -  tzokR  •  v 
3 


d\> 

dt 


The  last  term  in  Eq.  (35)  is  rate  of  elastic  strain  energy,  T ,  Specific  rate  of  strain  energy  is  in  the 
same  form  of  Eq,  (9),  except  the  effective  stress  for  linear  clastic  response  is 


CT  =  CT  y 


for  s  <  £r . 


(41) 


because  incompressibility  is  retained  in  considering  elastic  strain  energy,  strain  energy  rate  follows 
that  of  Eqs,  (13)  and  (14).  That  is 


4'  =  f  <jy  f{s  p )(— ) : 2  —  dV  =  ZttR  Vrv  r  °r , 

&  r  r  £  r 


t‘  £  dr 


(42) 


Again,  incompressibility  allows  us  to  write  —  =  — \ 

v  2 


ds 


.  Since  £  (<  ey)  is  a  very  small 


exp(^)-l 


number;  exp(-£) « 1  +— £  .  Therefore,  —  =  — .  Replacing  this  expression  in  Eq.  (42), 

2  2  r  3  £ 


'lJ=-&xR 


7jR2(jyV  I  where  the  proper  integration  limits  are  used.  Therefore, 

*V  3 £y  3 


•  8  , 

Y  =  -;c R'VyV 


(43) 


Eqs.  (36),  (40)  and  (43)  give  the  formula  of  each  term  involving  local  energy  in  elastic  zone. 
Placing  them  in  Eq.  (35),  we  obtain  the  radial  stress  on  yield  surface: 


<7 


D 


,  1  _  dv  1  5 , 

PXR—  +  -  pxv  ) 
3  dt  2* 


(44) 
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Combining  Eq.  (44)  with  Eq,  (IS),  we  obtain  the  governing  equation  for  dynamic,  transient  void 
growth.  The  development  of  a  dynamic  model  of  spherical  void  growth  is  thus  completed. 


3.  Threshold  condition  and  quasistatk  void  growth 

We  first  examine  the  situations  that  dynamics  effects  are  negligible.  Many  prior  studies  have 
been  devoted  to  examine  this  subject,  see  review  in  Section  l .  Among  others,  quasistatic  growth  of 
spherical  void  has  been  studied  by  Huo  et  al  (1999)  who  considered  surface  energy,  Wu  et  al  (2003) 
also  worked  on  the  model  with  materials  possessing  various  strain  and  strain- rate  sensitivity. 
However,  that  work  does  not  consider  the  effects  of  surface  energy.  The  example  calculations 
presented  in  the  following  only  involve  non -straining  materials;  applications  of  the  developed 
model  to  more  general  material  behaviors  will  be  reported  elsewhere. 

By  Eqs.  (33)  or  (44),  wc  recognize  that  radial  stress  on  yielding  surface  in  quasi  static  growth 


is 


(45) 


Eq.  (18)  is  then  simplified  to  the  following  equation  if  dynamic  effects  are  dropped  altogether: 


(46) 


Eq.  (46)  is  in  the  same  form  as  that  of  Huo  et  al  (1999)  except  the  lower  limit  of  integration  is  sY 
instead  of  zero  in  Hou  et  al  (1999),  The  lower  limit  of  zr  is  the  strain  al  boundary  of  plastic  zone  (r 
=  D),  Integration  of  Eq,  (46)  is  carried  out  for  non-strain  hardening  materials;  i.e.,  f{e)  =  1  -  By- 
changing  variables,  the  analytical  form  of  integration  is: 
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ds 


(47) 


exp<-s)-l 


=  |lnp-exp{-|s)], 

3  L 


By  integrating  between  aforementioned  limits,  we  obtain  the  relation  between  applied  pressure,  Pt 
and  void  size,  R: 


D  2y  2  2 

P-  —3 — <7y+  —  &r  In 

R  3  1  3 


i-4 

R 


I-exp(-'%r) 


(48) 


It  is  surprising  that,  after  many  papers  published  in  literature,  the  dosed  form  of  hq.  (48)  is  obtained 

3  3 

for  the  first  time.  By  approximating  1  -  exp(-  -  sY )  ~  —  sr  , 


P_ 

<*r 


U+!+^[o4)-1. 

R  aT  3  3  R  3ey 


(49) 


Obviously,  the  denominator  becomes  zero  iflowcr  limit  of  integration,  er ,  is  zero.  Thus,  the  model 
of  Huo  ct  al.  (1999)  breaks  down  for  single  void  embedded  in  infinite  matrix  (at  zero  void  fraction, 
Yf—  0). 


We  first  consider  a  hypothetical  case  in  which  surface  energy  can  be  ignored  ( y  -  0 ).  For 
instance,  systems  with  large  isolated  voids  fall  into  this  category.  Since  the  logarithmic  term  on  the 
R.HS  of  Eq,  (49)  represents  energy  dissipation,  it  cannot  be  negative.  At  initial  yielding,  zero 


R 3  2 

dissipation  implies  ln[(l  — - ]  =  0 :  and  the  corresponding  macroscopic  stress  by  Hq.  (49)  is 

R' 

P  -  2uy  /3  -  That  is,  void  starts  expansion  by  yielding  on  void  surface  at  P  =  2 vr  /3 .  The  same 
conclusion  can  also  be  obtained  by  examining.  Eq.  (45).  On  void  surface  where  initial  yielding 
starts,  the  radial  stress  is  zero,  cr.  =  <rD  =  0 .  By  Eq.  (45),  the  corresponding  macroscopic  stress  is 
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P  =  2o y  /3.  At  the  minimum  macroscopic  stress  of  P  —  2ar  /3for  initial  yielding,  the  condition  of 
R*  2 

3n[(l - r  ) - ]  =  0  gives  the  new  void  size  upon  initial  yielding  as, 

R  afj.- 


A  =  (l-^L)-ins|  +  fL 
R„  2  2 


(50) 


At  a  typical  value  of  £y  =  0,002 ,  the  growth  of  void  size  is  about  0,1%,  Although  0.1%  is  a  small 
expansion,  repeated  loading  can  raise  the  size  indefinitely.  For  example,  1 ,000  cycles  of  loading  at 
minimum  stress  of  P  =  2 <rr  /3  result  in  100%  void  expansion.  On  the  other  hand,  at  this  absolute 
minimum  loading,  the  growth  rate  is  zero  (v  =  0)  and  void  does  not  grow  at  all  in  theory.  At  a 
loading  below  the  minimum  stress,  no  plastic  yielding  ever  takes  place,  and  the  system  remains 
entirely  elastic. 


For  the  hypothetical  case  of  zero  surface  energy,  Fig.  2  shows  the  calculated  results  of 
P/(7y  versus  normalized  void  radius,  Rf  Rc ,  according  to  Eq.  (49)  for  a  range  of  yield  plastic 
strain  £y .  The  results  in  Fig.  2  are  supposed  to  be  identical  to  that  of  Huang  et  ah  (1991  )r  As  shown, 
equilibrium  void  size  increases  very  rapidly  as  the  applied  pressure  approaching  an  asymptotic 
value.  Equilibrium  void  size  does  not  exist  for  stress  exceeding  a  so-called  threshold  stress,  P 
(e.g„  Huang  et  al.,  1991 ;  Wu  et  ah,  2003),  The  threshold  stress  is  obtained  by  letting  R  &  in  Eq, 
(49);  i.e., 


p»=p | 


2  2 

=— (J  +  ln - ){Ty  . 

3  3 By 


(51) 


The  threshold  stress  of  Eq.  (51)  is  identical  to  that  of  Huang  et  ah  (1991)  and  others.  Therefore,  for 
a  loading  P  <  ,  void  grows  to  a  finite,  equilibrium  size  determined  by  Eq.  (49)  (shown  in  Fig.  2). 

Once  loading  exceeds  the  threshold  stress  P. ,  void  can  grow  indefinitely  with  growth  rate 
determined  by  the  inertial  associated  with  dynamic  growth  (see  Section  4). 
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With  zero  surface  energy  (y  =  0 )  in  Eq.  (49),  the  conditions  R  — >  ac  and  Ro  — ►  0  lead  to  the 

same  results  of  Eq.  (51).  The  latter  is  called  cavitation  limit  (c.g.,  Huang  et  aL,  1991).  However,  in 
light  of  the  complete  formulation  of  Eq.  (49)  with  effects  of  surface  energy,  the  conditions  R  ae 

2  y 

and  R.,  — >  0  are  not  the  same.  In  the  latter  case,  surface  energy  term, - ,  becomes  infinity  for 

R  cy  ' 

void  to  grow  from  R„  0  .  It  is  not  possible  for  cavitation  to  take  place  from  a  mathematical  point 
of  zero  size.  This  conclusion  is  quite  different  from  previous  studies  in  this  respect. 

Void  growth  with  effects  of  surface  energy  is  shown  in  Fig,  3,  For  comparison  purposes, 
material  data  ofHuo  ct  al  (1999)  arc  used  here;  i.e.,  y  =  1  J/nr,  <rr  =  200  MPa,  E  =  130  GPa,  and 
mass  density  p  =  8,950  Kg/nr.  Since  surface  energy  is  significant  only  for  very  small  voids,  initial 
void  radius,  R 0 ,  as  small  as  2  nm  are  shown  in  the  figure  (although  the  applicability  of  continuum 

approach  becomes  highly  questionable  at  this  small  length  scale).  Overall,  surface  energy  adds 
additional  requirement  on  applied  stress  to  expand  voids.  As  void  grows  much  bigger,  the  term 

2  v 

involving  surface  energy,  -  — ,  in  Eq.  (49)  diminishes;  the  ratio  PZ  o\.  approaches  the  asymptotic 

R  Uy 

value  of  Eq.  (51), 

Unlike  the  cases  with  y  -  0  which  is  monotonic,  an  interesting  characteristic  in  the  curves 
of  PZ  oY  versus  R  Z  Ro  in  Fig.  3  is  the  existence  of  peaks.  For  example,  for  Re  =  5nm  ,  Pic ty 
reaches  the  peak  value  of  5,838  at  RZ  R0  =1.328,  although  it  takes  lower  stress  P  for  the  void  to 
continue  growing  once  void  size  pass  the  peak.  Therefore,  if  the  applied  stress  exceeds  the  peak 
stress  of  P  /  <ty  =5.838,  void  becomes  unstable  and  expansion  continues  indefinitely.  We  will  call 
this  as  the  critical  stress ,  ex.. ,  for  unbounded  void  growth.  The  existence  of  critical  stress  is  not 

unexpected.  Void  growing  from  small  initial  size  needs  to  overcome  the  resistance  due  to  surface 
energy  (surface  tension),  which  decreases  with  large  void  size.  On  the  other  hand,  resistance  due  to 
clastic-plastic  deformation  increases  with  void  size.  Therefore,  there  is  a  transitional  void  size  at 
which  applied  stress  is  balanced  by  surface  energy  and  plastic  deformation,  At  that  critical 
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condition,  void  docs  not  grow.  The  critical  condition  is  corresponding  to  dP  /  dR  =  0  ,  where  P  is 
given  by  Eq.  (49).  It  can  be  shown  that  the  critical  condition  leads  to  the  following  equation: 


— tfr=Ar-i. 

y  R 


where  R  is  void  size  at  critical  condition.  Once  the  critical  void  size,  Rc  ,  is  obtained,  the 
corresponding  critical  stress,  o\. ,  is  calculated  by  Eq,  (49).  The  above  equation  is  solved 
numerically;  and  results  are  shown  in  Fig.  4  for  a  range  of  initial  void  radius,  Ra .  The  influence  ot 
specific  surface  energy  y  is  also  shown.  For  the  void  of  initial  size  Rf.  =  5nm  embedded  in  model 
material  with  <ry  =  200  MPa,  the  critical  stress  is  1.16  GPa  (5,838  ar )  for  y  =  1 ,  and  1,47  GPa 
(7.349  (7y )  for  y  =  2,  Loading  at  the  cavitation  limit  stress  of  Eq,  (51 )  does  not  lead  to  cavitation  at 

all,  due  to  the  rapidly  increasing  surface  energy  term  — — .  It  does  not.  lead  to  unbounded  growth 

R 

from  a  small  initial  void  either.  For  a  void  to  grow  indefinitely,  the  applied  stress  needs  to  overcome 
the  aforementioned  critical  stress,  er  ,  It  is  also  noted  that  the  critical  stresses  corresponding  to 
different  initial  void  sizes  eventually  converge  to  a  single  asymptotic  value  at  large  void  size,  which 
is  the  threshold  stress  of  Eq.  (51), 


With  the  effects  of  surface  energy  being  considered,  initial  yielding  takes  place  at  a  stress 
higher  than  P  =  2crr  /3  quoted  earlier.  According  to  Eq.  (49),  the  minimum  stress  for  initial  yielding 
would  be 


p=2y+2a„.,  (52) 

R  3 


Again,  the  same  expression  of  Eq.  (52)  can  be  obtained  by  examining  Eq,  (45).  On  void  surface 

2y 

where  initial  yielding  takes  place,  surface  energy  leads  to  a  surface  normal  traction  of  T  -  — , 

R 

pointing  toward  void  center  (Gent  and  Tompkins,  1969).  That  is,  the  radial  stress  on  void  surface  is 
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(j r  =  crn  =2 y/R.  By  Eq,  (45),  the  corresponding  loading  is  P=  2y  fR  +  2(Jr  / 3  .  This  expression  is 
identical  to  Eq,  (52).  The  dependency  of  macroscopic  stress  for  initial  yielding  on  void  size  is 
calculated,  and  sho  wn  in  Fig.  4  using  data  of  model  material.  Obviously,  the  effects  of  surface 
energy  are  dominant  tor  the  range  of  void  size  considered  in  the  figure.  For  void  of  radius  R  =  30 
ran,  contributions  of  surface  energy  and  yield  stress  to  the  loading  for  initial  yielding  are  equal  (for 
the  model  material);  at  larger  void  size,  the  role  of  surface  energy  becomes  less  significant, 

4.  Dynamic  void  growth 

Dynamic  void  growth  has  been  the  subject  of  a  number  of  papers,  including  Carroll  (1972), 
Cories  el  al  (1992,  1995),  Ortiz  and  Molinari  (1992),  Wu  et  al.  (2003).  The  present  model,  as 
represented  by  Eq,  (18)  in  conjunction  with  Eq.  (33)  or  (44),  however,  contains  correction  factor  a 
to  account  for  the  limited  plastic  zone  surrounding  a  growing  void;  and  fc  for  compressible  elastic 
zone.  In  light  of  these  improvements,  the  results  are  expected  to  be  more  accurate  than  those  in 
existing  literatures.  Again,  this  study  only  examines  non-strain  hardening  materials:  more  general 
applications  of  the  developed  model  will  be  reported  elsewhere.  The  considered  loading  is  a 
constant,  hydrostatic  stress  that  is  applied  instantly  at  time  zero.  This  type  of  loading  prevents  the 
additional  effects  of  loading  speed  like  that  considered  in  Wu  et  al.  (2003).  In  the  following,  steady- 
state  growth  will  be  examined  first,  followed  by  a  brief  discussion,  and  example  calculations  of 
transient  growth. 

4.1  Steady-state  growth 

In  steady  state  growth,  dRfdi  =  v  —  constant.  The  two  approaches  of  Sections  2.2  and  2.3 
lead  to  different  forms  of  radial  stress  on  yield  surface,  aD  in  Eq.  (33)  and  (44),  respectively. 

Before  engaging  a  full-flange  comparison  between  Eqs.  (33)  and  (44),  it  will  prove  rewarding  to 
examine  the  asymptotic  condition  of  void  growth  at  very  low  rate  (v  — >  0)  or  in  materials  with  near- 
incompressibility  I  c.  =*ao ).  In  either  ease,  we  have  the  condition  of  v  =  vfce  — >  0.  Replacing  aD 
from  Eq.  (33)  into  Eq,  ( ]  %}  for  steady  state  growth,  and  assuming  J?  >> 
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(53) 


^(a+my  =  P-P„, 

where  threshold  stress  P.  is  given  by  Eq.  (5 1).  1 1  can  be  shown  by  using  Eq.  (34)  that 

|  4  |  4 .y 

i v  —> - ( - )"  as  v  ->  0  ,  For  near-incompressible  materials  ( v  — >  1/21  — >  4/6  .  On  the  other 

ft  9  1-v 

hand,  replacing  &D  from  Fq.  (44)  into  Eq.  (18)  for  steady  state  growth: 


|(«  +  *:)v1  =  P-Pal 


(54) 


By  using  Eq.  (39),  it  is  easy  to  show  that  k  — >  3 f  0  as  v  ->  0 .  For  incompressible  materials,  the 


4  3 

coefficient  in  parentheses  on  the  LHS  of  Eq.  (.53)  is  cr  +  —  in  contrast  to  a  +  —  in  Hq.  (54). 


Accordingly,  the  model  of  Hunter  and  Crozicr  (1968)  gives  Lower  steady-state  growth  rate  in  near- 
incompressible  materials  than  that  predicted  by  the  model  of  Eq.  (54). 


For  near-incompressible  materials,  deformation  (clastic)  zone  extends  to  infinity  in  the  same 
way  as  has  been  addressed  in  the  literatures.  By  integrating  the  terms  on  the  RHS  of  Eq.  (2)  from  R 

to  infinity,  the  result  is  Actually,  the  governing  equation  for  near-incompressible  materials 

should  have  been  written  as: 


a 

*  ■  ■  V 

2 


P-4. 


(55) 


This  is  the  conventional  model  (e.g.,  Wu  et  al,,  2003,  with  added  surface  energy  term),  assuming 
incompressibility  for  material  surrounding  an  expanding  void.  If  one  replaces  the  relation  of 
a  =  3(1-1//?)  on  the  LHS  of  Eq.  (54),  the  equation  is  then  identical  to  Eq.  (55).  That  is,  the  model 
of  Section  2.3  reconciles  with  the  existing  dynamic  model  for  steady -state  void  growth  in  near- 
incompressible  materials.  Conversely,  the  model  of  Hunter  and  Crozicr  (1968)  fails  to  do  the  same. 
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To  further  compare  the  two  approaches  of  Sections  2,2  (i.e.,  Eq.  33)  and  2.3  {he.,  Eq.  44), 
we  are  focused  on  the  terms  co  in  Eq.  (34)  against  k  in  Eq.  (39).  The  calculated  results  of  these 
terms  (  {3k  and  (im )  are  shown  in  .Fig.  5{a)  for  two  values  of  Poisson’s  ratio  (0.25  and  0.499999) 
over  a  range  of  growth  rate  v.  The  nonnahzcd  growth  rate  v  =  0, 144  for  v  =  600  m/s  based  on  a 
Poisson's  ratio  v  =  0.25  .  As  growth  rate  v  approaching  zero,  the  value  of  /?*  approaches  to  3, 
regardless  of  the  value  of  Poisson^s  ratio,  v .  The  value  of  fko  does  not  converge  at  all  .  Since  the 
value  of  pK  deviates  from  3  at  higher  growth  rate  v,  the  difference  increases  between  the  present 
model  represented  by  Eq.  (54)  and  the  conventional  model,  of  Eq,  (55). 

Fig.  5(b)  shows  the  calculated  normalized  pressure,  (P-R.)/  pc] T  versus  growth  rate.  v. 
Material  data  are  taken,  from  the  model  material  except  the  indicated  Poisson's  ratio.  In  Fig.  5(b), 
the  curve  of  v  — >  0.5  by  using  the  present  model  of  Section  2.3  (labeled  Tszcng)  is  identical  to  the 
conventional  model;  i.c.,  Eq.  (55).  Incompressible  material  requires  higher  pressure  to  reach  the 
same  growth  rate;  due  to  larger  elastic  deformation  zone  (actually  infinity).  Although  the  results  of 
Tszeng  differ  only  slightly  from  that  of  Hunter- Crozier,  the  former  is  considered  more  favorable. 
The  seemly  big  discrepancy  between  k  and  &  in  Fig.  5(a)  docs  not  translate  to  significant 
difference  in  calculated  pressure  shown  in  Fig.  5(b),  The  reason  is  that  the  factor  a  is  much  larger 
than  jc  in  Eq.  (54)  or  O)  in  Eq,  (53);  changes  in  factors  k  or  to  do  not  significantly  affect  the 
pressure  P. 

According  to  Eq.  (54),  void  growth  rate  is  mainly  affected  by  the  threshold  stress,  ,  and 
less  by  the  relative  size  of  plastic  zone,  a  ,  That  is  the  case  when  growth  rate  is  modest.  Note  that 
the  threshold  stress,  ,  is  directly  related  to  yield  stress  via  Eq,  (51).  In  the  situation  of  P  »  Pif_, , 
however,  the  role  of  relative  size  of  plastic  zone,  a  ,  becomes  more  significant.  At  certain  point,  the 
higher  yield  stress  leads  to  a  smaller  value  of  a  so  much  so  that  the  growth  rate  becomes  higher. 

To  make  this  point.  Fig.  6(a)  shows  the  macroscopic  pressure,  P,  versus  void  growth  rate,  v,  for  the 
model  material  with  two  levels  of  yield  stress:  a,.  «=  100  MPa  and  200  MPa.  In  the  regime  of  low 
pressure  (therefore  low  growth  rate),  material  with  higher  yield  stress  grows  slower  for  the  same 
pressure.  At  a  pressure  P  greater  than  a  crossover  of  12.16  GPa,  the  contrary  is  true  (although  hard 
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to  read  from  the  figure).  That  is,  material  with  higher  yield  stress  grows  faster  for  the  same  pressure. 
By  using  Eq.  (54),  the  crossover  growth  rate,  v*,  can  be  determined  by 


v*  =  [- 


2 


■rfj] 


-r. 


(or,  +  a:  -a.-  - k, ) 


(56) 


In  this  expression,  subscripts  4TT  and  4TT  denote  the  reference  and  comparison  materials, 
respectively.  Since  the  factor  k  is  a  function  of  growth  rate,  Eq,  (56)  is  solved  numerically  for 
specified  material  data.  Results  of  crossover  growth  rate,  v*,  and  pressure,  P*,  are  shown  in  Fig. 
6(b)  using  our  model  material  as  the  reference  material  (crr  =  200  MPa);  comparison  material  has 
yield  stress  ranging  between  25  and  200  MPa.  At  a  pressure  P  <  P  \  (or  v  <  v*),  growth  rate  is 
lower  for  materials  with  high  yield  stress.  Conversely,  at  a  pressure  P  >  P*  (or  v  >  v*),  growth  rate 
is  higher  for  materials  with  high  yield  stress.  Since  Eq.  (56)  is  fully  symmetric  with  respect  to 
reference  and  comparison  materials,  the  roles  of  these  materials  can  be  swapped,  yet  crossover 
growth  rate  remains  the  same. 

At  this  point,  it  calls  for  the  remark  that  some  of  the  calculations  in  the  regime  of  high 
growth  rate  shown  in  Fig,  6(a)  and  (b)  are  questionable.  According  to  Eq.(37),  elastic  zone  is 
bounded  between  r  =  > D  and  cj .  At  a  growth  rate  v  -  clv  f  ft  t  yield  surface  reaches  the  dilatational 

wave  front.  According  to  Eq.  (39),  &  =  0  for  v  >  cc  i  ft  .  The  limiting  growth  rate  of  ce  /  ft  is  shown 

in  Fig.  6(b)  for  the  considered  range  of  yield  strength.  As  shown,  majority  of  the  calculations  fall 
beyond  the  limit.  In  order  to  better  handle  void  growth  at  higher  rate  or  high  loading,  the  model  has 
to  include  compressibility  in  plastic  zone,  and/or  allow  discontinuity  due  to  shockwave  (Hunter  and 
Crozier,  1968).  Since  k  =  0  for  v  >  ce  I  ft ,  it  is  interesting  to  observe  that  radial  stress  on  yield 
surface,  oy,,  is  the  same  as  that  in  quasistatic  growth  of  Eq.  (45). 


4,2  Transient  growth 


In  the  following  development,  the  approach  of  Section  2.3  is  employed.  The  complete 
governing  equation  is  obtained  by  Eq.  (IS)  with  aD  from  Eq.  (44): 


-p(a+K)R 
3  at 


(^r)2  =  P-  \\<ty+\er  ln[(l  -  4) 


d'~R  p(a  +  K ).dR^  _  (2  2  .  ...  Ri.  2 


at 


]+4r 


v- 


R  2s.,  R 


(57) 


This  equation  is  solved  numerically,  because  kt  is  a  function  of  average  growth  rate  too.  Only  brief 
applications  of  the  model  are  given  below,  further  discussion  will  be  reported  elsewhere. 


According  to  Section  3,  a  critical  stress  tr.  is  identified,  which  corresponds  to  the  peak  in 

the  curve  of  stress  versus  void  size  (see  Fig.  4),  Theoretically,  the  critical  stress  approaches  the 
threshold  stress  (or  so-called  cavitation  limit  stress)  of  Eq.  (51)  for  infinite  void  size.  Two  loading 
levels  are  considered  in  the  following  calculations:  subcritieal  (P  <  <JC )  and  supercritical  loading  (P 

>  <jc ).  The  critical  stress  is  based  on  the  initial  void  size,  like  that  shown  in  Fig.  4.  It  needs  to 

clarify  that,  the  work  of  Wu  et  ah  (2003)  also  considers  supercritical  and  subcritieal  loadings,  but 
their  critical  stress  is  threshold  stress  (or  cavitation  limit  stress)  m  our  terminology. 

For  subcritieal  loading,  void  starts  growing  once  the  RIJ.S  of  Eq.  (57)  becomes  positive.  It 
only  grows  to  a  limited  size  smaller  than  the  critical  size.  As  an  example,  for  Ro  =  5  nm,  the  critical 

stress  <JL  is  5.834f7p  and  critical  void  radius  R.  =  1.325/?.,  -  6.625  nm.  Fig.  7(a)  shows  the  growth 
of  the  void  subjected  to  a  subcritieal  loading  P  =  5J<rr  <5.834 crY  .  The  growth  rate  increases  at  the 
moment  load  is  applied:  it  reaches  a  maximum  shortly,  and  starts  declining  to  zero.  The  void  ceases 
growing  at  R  =  6.5  nm,  which  is  less  than  the  critical  radius  Rc  of  6.625  nm. 


Fig.  7(b)  shows  the  results  when  a  supercritical  loading  P  =  6a  Y  >  5,83  8cr,,  is  applied.  Void 
continues  growing  beyond  the  critical  size  which  is  6.625  nm;  and  is  expected  to  grow  unbounded  if 
the  constant  loading  maintains  unchanged.  The  grow  rate  is  not  monotonic  though;  it  reaches  a 
plateau  before  resumes  increasing  trend  in  the  small  time  frame  of  0,1  ns.  The  temporary'  reduction 
in  growth  rate  is  due  to  the  growing  resistance  of"  elastic-plastic  deformation  (which  is  zero  in  the 
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beginning  of  growth).  The  critical  condition  is  marked  in  Fig.  7(b),  which  indicates  that  largest 
resistance  (relative  sense)  against  void  growth  appears  at  critical  condition.  Fig.  8(a)  shows  growth 
rate  on  void  of  initial  size  tf.  =  5nm  subjected  to  a  range  of  macroscopic  stress.  As  shown, 
macroscopic  stress  higher  than  5.8<xr  leads  to  unbounded  growth  whereas  a  stress  of  5. Toy  does 
not.  However,  it  is  noticed  that  the  stress  of  5, Sen,  is  actually  lower  than  the  critical  stress  of 
5.834cy  corresponding  to  R0  -  5nm .  The  reason  unbounded  growth  is  found  at  a  stress  of 
5. Soy  <  o’,  is  that,  in  dynamic  growth,  void  size  may  exceed  the  critical  size  corresponding  to 
current  void  size  before  growth  ceases,  Growth  rate  for  voids  of  a  range  of  initial  size  is  shown  in 
Fig.  8(b).  The  loading  P  is  slightly  higher  than  the  critical  stress  cr.  corresponding  to  void  of 

R  =  5nm ,  but  sufficiently  greater  than  those  corresponding  to  Ay  =  Inm  or  higher  so  that  the 

aforementioned  plateau  docs  not  appear.  For  smaller  void,  the  growth  ceases  shortly  after  expansion 
starts.  The  numerical  results  of  Fig.  8  indicate  that,  for  supercritical  loading,  void  growth  rate 
reaches  a  steady  value  at  relatively  large  void  size,  This  is  discussed  in  previous  section, 

To  gam  a  better  understanding  of  the  void  growth  process,  energy  components  of  Eqs.  (5), 
(6).  (8),  and  (15)  for  the  plastic  zone  and  Eqs,  (36),  (40)  and  (43)  for  elastic  zone  are  calculated  by 
using  the  solved  growth  rate;  results  are  shown  in  the  form  of  area  chart  in  Fig.  9.  Each  of  the 
components  is  normalized  to  the  total  energy  supplied  by  the  remotely  applied  pressure.  Numerical 
results  indicated  the  sum  of  calculated  energy  components  by  Eqs.  (6),  (8),  and  (1 5)  is  equal  to  the 
total  energy  of  Eq.  (5)  for  plastic  zone'  similarly  for  elastic  zone.  This  confirms  that  the  equations 
are  solved  properly.  According  to  Fiq.  (8),  kinetic  energy  in  plastic  zone  has  two  components: 

2 npaRrv  for  inertia  and  (4/3 )npczR  'vdvUh  for  acceleration.  As  shown  in  Fig.  9(a)  for  a  pressure 
slightly  above  the  critical  stress  of  oy  =  5.1 5GPa  for  void  size  R,,  - 10 nm ,  surface  energy  and 

kinetic  energy  for  acceleration  represent  the  major  energy  expenditure  at  the  beginning.  Energy  in 
elastic  zone  is  not  negligible  either.  As  void  grows  bigger.,  surface  energy  and,  particularly,  kinetic 
energy  associated  with  acceleration  drop  quickly  At  the  same  time,  plastic  dissipation  and  kinetic 
energy  associated  with  inertia  in  plastic  zone  increase  rapidly.  Energy  in  elastic  zone  mainly  spent 
on  clastic  strain  energy  which  maintains  constant  proportion;  kinetic  energy  represents  a  very  small 
portion.  As  void  becomes  even  bigger,  energy  expenditures  are  divided  among  plastic  dissipation, 
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kinetic  energy  associated  with  inertia  in  plastic  zone,  and  elastic  strain  energy.  At  higher  loading,  P, 
Fig.  9(b)  shows  higher  proportion  of  kinetic  energy  associated  with  inertia  in  plastic  zone,  and 
clastic  zone  as  well 

The  proportion  of  energy  expenditure  in  plastic  and  elastic  zones  maintains  constant  as  void 
becomes  bigger.  According  to  Eq.  (5),  ratio  of  energy  input  to  plastic  zone  to  total  energy  provided 
by  macroscopic  stress  P  is  <xD  IP.  As  void  grows  under  constant  macroscopic  stress  P,  growth  rate 
reaches  a  steady  state  value.  The  radial  stress  on  yield  surface,  ,  becomes  a  constant,  according 
to  Eq.  (44).  This  explains  the  observed  constant  proportion  of  energy  in  plastic  and  elastic  zones. 

We  can  apply  the  similar  treatment  to  all  other  energy  components  as  well. 

5.  Cond tiding  Remarks 

A  mathematical  model  is  established  to  properly  account  for  material  compressibility  and 
limited  size  of  elastic  deformation  in  dynamic  growth  of  spherical  void.  It  corrects  the  issue  of 
excessive  kinetic  energy  expenditure  in  elastic  zone.  A  list  of  most  signi  ficant  remarks  emerging 
from  the  study  is: 

1 .  The  steady- slate  solution  of  radial  stress  on  yield  surface  obtained  by  Hunter  and  Crozier  (1968) 
does  not  reconcile  with  the  solution  assuming  incompressibility.  On  the  other  hand,  an  approximate 
solution  is  developed,  which  can  be  used  for  general,  transient  growth. 

2.  A  closed  form  of  the  stress -void  size  relation  in  quasi  static  growth  is  obtained  for  the  first  time.  It 
shows  initial  yielding  on  void  surface  lakes  place  at  P  —  2/3  Y  for  large  voids.  For  smaller  voids  of  a 
few  nanometers  in  size,  macroscopic  stress  needs  to  overcome  a  critical  stress  for  unbounded 
growth.  Critical  stress  results  from  higher  surface  energy  at  decreasing  void  size.  Based  on  the  same 
reason,  cavitation  from  void  of  zero  size  does  not  seem  to  be  possible. 

3.  A  crossover  condition  is  identified  pertaining  to  the  effects  of  yield  stress  on  void  growth  rate.  At 
a  growth  rate/stress  below  the  crossover  condition,  higher  yield  stress  (i.o,.  work  hardened)  leads  to 
lower  growth  rate.  Reversal  is  true  at  growth  rate/stress  exceeding  the  crossover. 
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Figure(s) 


Figures 


Fig,  h  Schematics  of  a  expanding  spherical  void  embedded  in  infinite  matrix,  subjected  to 
remotely  applied  stress  P.  Void  radius  i?,  plastic  zone  radius  D.  Dilatational  elastic  wave  propagates 
at  speed  clv.  Initial  void  radius  is  R0 . 


Fig.  2.  Macroscopic  pressure  versus  void  size  in  quasistatic  growth.  sr  -  aY  i  E  is  total  strain  at 
yielding.  Surface  energy  density  y  is  zero. 
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Fig.  3.  Macroscopic  pressure  versus  void  size  in  quasistatic  growth  of  void  from  a  range  of  initial 
size..  Ro .  Surface  energy  density  y  -  L 


(a) 


(b) 


Fig,  4,  (a)  Critical  stress,  and  (b)  critical  void  radius  for  unbounded  void  growth  as  functions  of 
initial  void  size.  The  critical  stress,  is  normalized  by  the  yield  strength,  <jy  =  200 MPa  for  the  model 
material.  Also  shown  in  (a)  are  the  pressure  for  initial  yielding.  Two  levels  of  surface  energy 
density  are  considered, 
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(a) 


Fig.  5.  (a)  Factors  k  and  0)  as  functions  of  steady-state  void  growth  rate;  (b)  Normalized  pressure,, 
(P— p  )/  pc l ,  as  function  of  steady-state  growth  rate,  v.  tor  model  of  Section  2.2  (Hunter/Crozicr) 
and  Section  2.3  (Tszeng).  Property  data  are  taken  from  model  material,  except  twro  indicated  values 
of  Poisson's  ratio  arc  used,  v  -  0,499999  for  near-incompressible  material. 
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Fig,  6,  (a)  Pressure,  P,  as  function  of  steady-state  growth  rate,  v .  for  modd  materials  but  two  values 
of  Young’s  modulus  as  indicated.  Crossover  is  found  at  v*=  1,000.16  m/s  and  ?*=  12.16  GPa.  (b) 
Crossovergrowth  rate,  v*,  and  pressure,  P*,  as  function  of  yield  stress,  <Tr ,  Model  material  (<Jr  — 
200  MPa)  is  used  as  reference.  Limiting  growth  rate  corresponding  to  a:  =  0  (zero  elastic  zone)  is 
also  shown. 
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Fig,  7,  Dynamic  growth  of  void,  (a)  Suhcritical  loading,  P  =  5.1ay  <5.834<ry  .  (b)  Supercritical 
loading,  P  -  ba,  >  5.834crr  .  Initial  void  radius  R0  =  5  nm.  Property  data  from  model  material. 
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(a) 
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Fig.  8,  Growth  race  in  subcritical  and  supercritical  loading,  (a)  Loading  stress  is  below  and  above 
the  critical  stress  of  er:.  =  5.834c7r  ,  Initial  void  size  Ra  =  5nm  ;  (b)  Growth  from  void  of  a  range  of 

initial  size.  P  —  6a Y  >  5,834crK . 
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Fig,  9 .  Area  chart  of  energy  expenditure  in  the  process  of  supercritical  void  growth.  All  energy 
components  are  normalized  to  the  total  energy  supplied  by  the  externally  applied  pressure.  Initial 
void  size  Rt,  —  10  nm.  (a)  P  =  6etf  ,  (b)  P  =  10tfF . 


